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Abstract 



We present a lattice QCD calculation of the scattering phase shift for the I = 
2 5- wave two-pion system using the finite size method proposed by Liischer. 
We work in the quenched approximation employing the standard plaquette 
action at /3 = 5.9 for gluons and the Wilson fermion action for quarks. The 
phase shift is extracted from the energy eigenvalues of the two-pion system, 
which are obtained by a diagonalization of the pion 4-point function evaluated 
for a set of relative spatial momenta. In order to change momentum of the 
two-pion system, calculations are carried out on 24^ x 60, 32^ x 60, and 48^ x 60 
lattices. The phase shift is successfully calculated over the momentum range 
< p2 < 0.3 GeV^. 



PACS number(s): 12.38.Gc, ll.15.Ha 
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I. INTRODUCTION 



Calculation of scattering phase shift is an important step for expanding our understand- 
ing of strong interactions based on lattice QCD beyond the hadron mass spectrum. For 
scattering lengths, which are the threshold values of phase shifts, several studies have al- 
ready been carried out. For the simplest case of the two-pion system, the 1 = 2 scattering 
length has been calculated in detail [1-7] including the continuum extrapolation [5-7]. There 
is also a pioneering attempt at the / = scattering length [2], which is much more difficult 
due to the presence of box and disconnected contributions. For the scattering phase shift, 
in contrast, there has been only one calculation for / = 2 by Fiebig et.ai, who used lattice 
simulations to estimate the effective two-pion potential and used it to calculate the phase 
shift in a quantum mechanical treatment [8]. 

In this article, we calculate the 1 = 2 S-wave two-pion scattering phase shift applying 
the Liischer's finite size method [9,10]. Technically the key feature is the extraction of the 
two-pion energy eigenvalues from the pion 4-point function. This is successfully solved by a 
diagonalization method proposed by Liischer and Wolff [11] for 0(3) non-linear a model in 
2-dimensions. We also extract the scattering length from the phase shift data, and compare 
it with previous calculations. We work in quenched lattice QCD employing the standard 
plaquette action for gluons and the Wilson fermion action for quarks. 

We wish to mention that the study of the two-pion scattering phase shift also has im- 
portant impact on the calculation of the K ^ im decay amplitudes. A direct calculation of 
the amplitude from the 4-point function {0\7i (1,^)71 (tT^)Hw{tH)K(tK)\0) is very difficult, as 
pointed out by Maiani and Testa [12], because the 4-point function at large times is domi- 
nated by the two-pion ground state with zero relative momenta, which differs from the final 
state of the decay having a non-zero relative momentum. An exception is the amplitude 
from the K meson to the two-pion ground state itself, because this can be calculated by 
taking the two-pion state with zero relative momentum in the final state. However, the 
amplitude thus obtained is unphysical, and a reconstruction of the physical amplitude using 
some effective theory of QCD, for example chiral perturbation theory (CHPT), is needed. 
Using such an effective theory causes large uncertainties in the lattice prediction of the decay 
amplitude. Hence, a method for direct calculation of the — > tttt decay amplitude has been 
strongly desired. 

Recently Lellouch and Liischer [13] obtained a relation between the lattice and the phys- 
ical amplitude in the two-pion center of mass system with the energy E^^ = nix- In their 
derivation no effective theory is used. Lin et.al. [14] derived the relation from a different 
approach, and extended it to the general two-pion system with the energy E,^ ^ ttik- They 
also investigated the limitation of the relation. 
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In order to apply the relation to obtain the physical decay amplitude, one has to calculate 
the amplitude from K meson to the two-pion energy eigenstate with non-zero momenta on 
the lattice. This is the same problem as one encounters in the calculation of phase shifts 
using the Liischer's method. Thus study of the two-pion system represents a first step toward 
K — > TTTT decay. 

This paper is organized as follows. In Sec. II we describe the formahsm for calculation 
of the scattering length and phase shift [9,10]. We also discuss the method of extraction of 
energy eigenvalues of the two-pion system from the pion 4-point functions. The simulation 
parameters used in this work are given in Sec. III. In Sec. IV we analyze the behavior of 
the 4-point functions, and show that the diagonalization technique proposed by Liischer and 
Wolff allows to extract the energy eigenvalues. We then present results for the pion phase 
shift. Our conclusions are given in Sec. V. A preliminary report of the present work was 
presented in Ref. [15]. 



II. METHODS 
A. Finite Size Method 



The energy eigenvalues of a non-interacting two-pion system on a finite periodic box of 
a size are quantized as follows : 

K = 2-^m2+p2 , p^ = (27r/Lf -n, neZ. (1) 

In the interacting case the n-th energy eigenvalue is given by 



E^ = 2-^ml+Pl, pl={2n/Ly-n, n^Z, neZ. (2) 

The energy eigenvalue is written as that of the non-interacting two-pion system with mo- 
mentum pn and — Pn, but the quantity n — L^/(27r)^ -p^ is not an integer. The momentum 
satisfies the Liischer relation [9,10], 

Un5{Pn) = , . , (3) 
where S{pn) is the 5"- wave scattering phase shift at infinite volume and 

Zoo(fc;n) = -^. (4) 

Using (3), we can obtain the scattering phase shift from the energy eigenvalue calculated in 
lattice simulations. The scattering length is given by gq = limp^o tan 6{p) /p. 
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In the limit of large volume or weak two-pion interactions, we find 



^2 -''-0{1/V') 



Pn-Pn 



or n — n = 0{1/L) 



(5) 



from (3) and (4). Therefore, taking the volume to be large in lattice calculations, we can 
employ an expansion of Zoo{l;n) around n e Z given by 



Att ■ Zoo(l;n) 



N 



+ ji"^ ^00 (j; n)-{n- n) 



n — n N-Kx 



(6) 



where 



Zoo{j;ri) = lim 



An ■ Zoo{j; n)- Nn-{n- n)~ 



(7) 



and = I^mgzs 5(m^ — n). In this work wc use the expansion (6) with N = 10, with which 
the numerical error for all our simulation parameters are under 0(10^^). The numerical 
calculation of ^oo(j ; ^) is discussed in Ref. [9] . The values for several j's and n's are tabulated 
in Table I. 



B. Extraction of Energy Eigenvalues of Two-Pion System 

In order to obtain the energy eigenvalues of the two-pion system we construct the pion 
4-point function : 

Gii^Kt) = m^m(^-\tsm . (s) 

Here fln{t) is an interpolating field for the 5"- wave two-pion system at time t given by 

^n{t) = ^■Y.ARiPn),tM-RiPn),t) , (9) 
R 

where 7i{pn,t) is the pion interpolating field with lattice momentum at time t. The 
vector Pn satisfies pf^ = (27r/L)^ ■ n {n E Z), and R is an element of cubic group which has 
48 elements. The summation over R is the projection to the sector of the cubic group, 
which equals the S'-wave state in the continuum ignoring effects from states with angular 
momentum L > 4. 

For source we use another operator 1^^^^-' (t) defined by 

1 Nr 

^^i^^H^) = TT ■i:<Pn,t,Q7ri-p^,t,r],) , (10) 

where 
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y 



(11) 



The field n{f)n, t, Vj) is defined as 7r(p„, t, ^j) by changing C,j{x) to rij{x). The functions ^^(a;) 
and r]j{x) are orthogonal complex random numbers in 3-dimensional space, whose property 
is 
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Nr 



ir ■ E = - y) ■ (12) 



The pion 2-point function is constructed as 

Nr 



1 i^R 

Gl^'^'^Kt) - ^ ■ E (0|7r(p„,t)7r(-p-;,t5,e.)|0) ■ (13) 



When the number of random noise source is taken large or the number of gauge config- 
urations becomes large, we expect 

^ Gl{t) = (0|7r(p-;,t)7r(-p-;,i5)|0) , (14) 

and the 4-point function will be symmetric under exchange of the sink and source momenta. 
In our numerical calculations we use U(l) random numbers and take = 2. The number 
of configurations is 200, 286 and 52 depending on the lattice size as shown in Sec. III. We 
always check the symmetry of the 4-point function across the midpoint in the temporal 
direction before analysis. 

The 4-point function can be rewritten in terms of the energy eigenvalue Ej and eigenstate 
\Clj) as 



(o|ajn,)(n, |Q^|o) 

where fi„ = ^^n(0) and we assume non-degeneracy of energy eigenstates. The j-th energy 



Ej — 2 - ■sjm'l + satisfies the Liischer relation (3). Since the matrix element {0\Qn\^m) is 
not diagonal generally, the 4-point function Gnm{t) contains many exponential terms and is 
not a diagonal matrix with respect to the momentum indices n and m. For simplicity we 
introduce the following matrices : 



Ar,m{t)-Snm-e-^-<'-''^ (16) 



and rewrite the 4-point function in the following matrix form. 
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G{t) = V A{t) , 



(17) 



where n and m are regarded as matrix indices. 

The extraction of the energy eigenvalues from muhi-exponential Green function such as 
(17) is non-trivial. One can attempt multi-exponential fitting to extract them, but it is very 
difficult in general. A method of extraction was proposed by Liischer and Wolff [11]. They 
applied it to the 0(3) non-linear a model in 2-dimensions and obtained the scattering phase 
shift. This method has been used for many statistical systems [16] and also for the 1 = 2 
two-pion system of QCD [8]. In their method the following matrix is diagonalized at each t, 



M{t,to) = G{to)-'/' G{t) G{to) 



-1/2 



(18) 



where to is some reference time. The eigenvalues X{t,to) of M{t,to) can be obtained easily 
from (17) and (18) by 



X{t,to) = Ey 



M{t,to) 



Ev 



G{t)G{to)-' 



Ev 



VA(t)A(to)-^V-^ 



exp ( —Ej ■ (t — t[ 



Ev 



J=0,l,2^ 



(19) 



Therefore after diagonahzation of M{t,to) we can obtain the energy eigenvalues Ej by a 
single exponential fitting. 

In actual calculations we can not calculate all the components of the 4-point function 
precisely. We have to set a momentum cut-off p^^^ = (27r/L)^ • N. Here we expect that the 
components of Gnm{t) for n,m <k are dominant for the k-ih. eigenvalue \k{t) in the large 
t and to region, while the components n,m > k are less important. In this work we set 
and t large and investigate the cut-off dependence ior N > k . 



III. SIMULATION PARAMETERS 



Our simulation is carried out in quenched lattice QCD employing the standard plaquette 
action for gluons at /3 = 5.9 and the Wilson action for quarks. Quark masses are chosen to 
be the same as in the previous study of the quenched hadron spectroscopy by CP-PACS [17], 
i.e., n = 0.1589, 0.1583, 0.1574, and 0.1566, which correspond to m^/mp = 0.491(2), 
0.593(1), 0.692(1), and 0.752(1). The lattice cut-off is estimated from the p meson mass, 
and equals l/a = 1.934(16)GeV. 

In order to examine finite-size effects for the scattering length and to change the momen- 
tum for the phase shift, lattice simulations are carried out for three lattice sizes with a fixed 
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temporal size T — 60. The number of configurations and the momentum = (27r/L)^ • n 
for each lattice size are tabulated below. 

configurations n 

24^ 200 _Q_ , i , 2 

32^ 286 _Q_ , i , J_ , 3 

48^ 52 _Q_ , , J_ , 3 . (20) 

Here we calculate the phase shift at the momenta marked by under-bar; those un-marked 
are used to examine the momentum cut-off effects. The momenta in units of GeV^ chosen 
in this work are plotted in Fig. 1. 

We note that the two-pion energy eigenstates are not degenerate for n < 6. Since the 
effects from the states n > 6 can be thought to be negligible for the first several low-energy 
states, the non-degeneracy assumption in the derivation of the diagonalization method in 
the previous section is justified. 

Gluon configurations are generated with the 5-hit heat-bath algorithm and the over- 
relaxation algorithm mixed in the ratio of 1 : 4. The combination is called a sweep and 
we skip 200 sweeps between measurements of physical quantities. Quark propagators are 
solved with the Dirichlet boundary condition imposed in the time direction and the source 
operator set at = 8 to avoid effects from the temporal boundary. 

IV. RESULTS 
A. Effects of Diagonalization 

In Fig. 2 we show examples of effective mass of the pion propagator G^{t) for momenta 
n = 0, 1,2 ( = (27r/L)^ ■ n ) at mj^/mp = 0.491 on a 32"^ lattice. The source operator 
is located at ts = 8. We observe a clear plateau over the time range t ~ 18 — 46 for small 
momenta, but the signal becomes noisier for large momenta. We also find very large effects 
from the temporally boundary for t > 46. 

The pion 4-point function Gnm{t) defined by (8) is plotted in Fig. 3 for the same param- 
eter. The signal is very clear, and we see that the off-diagonal elements (n 7^ m) are not 
negligible. This means that the overlap is not diagonal, i.e. Vnm 9^ ^nm in (15). Wc also 
observe that the 4-point function is almost symmetric under the exchange of the sink and 
source momenta, but the statistical errors are not symmetric. In the lower frame of Fig. 3, 
for example, Gi2{t) suffers from large statistical error, while that of 6*21 (^) is very small. In 
the following analysis we assume symmetry of the magnitude of error, and substitute the 
component with large statistical error by the symmetric partner with smaller error. We also 
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see evidence of the presence of many exponential terms in the lower frame of Fig. 3. The 
sign of Gi2{t) and G2i{t) is flipped at i ~ 36 — 38. This is possible only if more than two 
exponential terms are present. 

In order to examine the effects of diagonalization, we calculate two ratios defined by 

Rr,{t)=Gr,n{t)-[l/Gl{t)f , (21) 

Dn{t) = Xn{t, to) ■ [Gl{to)/Gl{t)] ' , (22) 

where A„(t,to) is the n-th eigenvalue of M{t,tQ) calculated with a finite momentum cut-off 
p^yt = {271 /Ly ■ N. If the 4-point function contains only a single exponential term, i.e. 
Gnm{t) oc 5nm ' ew[-En ■ {t - ts)]-, then 

Rn{t) = A ■ e-^^"-(*-*s) , (23) 

where Ai?„ = En — En and A is a constant. If the momentum cut-off is sufficiently large, 
then the eigenvalue behaves as A„(t, to) = exp[— A£^„ • {t — to)] and 

D^{t) = e-^^"-(*-*°) . (24) 

In these cases we can obtain the energy shift A£^„ = E^ — E^ easily from the ratio Rn{t) or 
Dnit) by a single exponential fit. 

In Fig. 4 the ratio Rn{t) and -D„(t) for the ground state n — are plotted for all quark 
masses and lattice sizes in this work. For Dn{t) the momentum cut-off pl^^ — (27r/L)^ • N 
is set at = 1 and the reference time is taken to be to — 18. We divide Dn{t) by a 
constant Dn{ts) to facilitate a comparison with Rn{t). The statistical errors are very small 
and the diagonalization does not affect the result. We also checked the momentum cut-off 
dependence by taking N = 2 and confirmed that it is negligible. In previous calculations of 
the scattering lengths [1-7] the ratio Ro{t) was used to extract the energy shift AEo- Our 
calculation demonstrates the reliability of these calculations. 

We compare the ratios for the first exited state n = 1 in Fig. 5. The momentum cut-off 
is set at N = 1 and N = 2. We divide Dn{t) by a constant Dn{ts) as for the case of n = 0. 
The diagonalization is effective for smaller quark masses and smaller lattice sizes, while it 
is less so for larger quark masses and larger volumes. The momentum cut-off dependence 
is negligible for all parameter region, however. We see a strange behavior near t = 36. We 
consider that this is either due to insufficient statistics or an effect of the temporal boundary. 
We then fit the ratio by a single exponential form over the time range consistent with the 
single exponential behavior. The fitting range for each parameter is tabulated in Table III. 

A similar comparison for n = 2 { the second exited state ) is made in Fig. 6. The mo- 
mentum cut-off is set at N — 2 and N — 3. We observe again that the diagonalization is 
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effective for smaller quark masses and smaller lattice sizes. The momentum cut-off depen- 
dence is small for all parameter region as for the case of n = 1. Compared with the n — 
and n — 1 cases, the signals are noisier. We observe a strange time dependence in the data 
at niT^/mp — 0.491 and 0.593 on a 32^ lattice at i ~ 30 — 46. For these data we restrict the 
fitting range to t = 18 — 32. We remove results at these parameters from our finial analysis. 
In other data clear signals of the single exponential behavior are seen for t > 18. The fitting 
range for each parameter is hsted in Table IV. 

Prom these results we conclude that the momentum cut-off should be taken N >nioic the 
energy shift A£'„. The results of the energy shift AEn obtained by the single exponential 
fitting of the ratio Dn{t) are tabulated in Tables II, III, and IV, where we take the 
momentum cut-off N — n, and the reference time to — 18. In the tables we also quote the 
scattering amplitude A{pn) defined by 

AW„) = . I , (25) 

where we normalize the amplitude as limp^o^(p) = " '^-k- 

B. Results of Scattering Length 

For n = the values of arc very small as shown in Table II. Therefore we may write 
yl(p„)/m^ ~ ao/'"^7ri and use results for n = to evaluate the scattering length. 

In Fig. 7 we recapitulate the recent results of JLQCD [6] and Liu et.al. [7] for the 1 = 2 
pion scattering length. The two values of Liu et.al. denoted as (Scheme I) and (Schema II) 
refer to their two different treatments of the finite volume corrections. The two values of 
JLQCD correspond to two different fitting functions for extraction of the energy shift from 
the ratio Ro{t), (LIN) used a linear fit in t while (EXP) employs a single exponential in t. 
Figure 7 shows that the lattice cut-off effect is strongly dependent on the choice of the fitting 
function. However, the dependence disappears toward the continuum limit. Compared with 
JLQCD the lattice cut-off effect of Liu et.al. is very small, since their calculation is carried 
out with an improved gauge and improved Wilson fermion action on anisotropic lattices, 
while the actions of JLQCD are the standard plaquette and the Wilson fermion actions. The 
values extrapolated to the continuum limit are consistent with the CHPT prediction [18] as 
shown in Table V. 

Since we use the same actions as those of JLQCD, we compare our results with theirs 
at the same gauge coupling constant /3 = 5.9 in Fig. 8. Here our data on a 48^ lattice are 
omitted, because those are consistent with the results on 24^ and 32^ lattices within very 
large statistical errors of those on the 48^ lattice (see Table II) . Our data for the scattering 
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length are different from those of JLQCD obtained by a linear fit (LIN) by about 2.5(7, 
whereas we find consistency among results obtained with the exponential fitting for four 
different lattice sizes, i.e. 24^, 32^, 48^ from the present work, and 16^ from JLQCD. In 
Fig. 8 we observe that both ours and JLQCD results at /3 = 5.9 are far from the CHPT 
prediction ao/rriT^ — —2.265(51) 1/GeV^. This is due to finite lattice cut-off effects, which 
are rather large for the standard actions as shown in Fig. 7. 

Here we comment on the choice of the fitting function for the ratio Ro{t). In our analysis 
we have assumed a single exponential behavior, i.e. Ro{t) ~ Z ■ exp{—AEo{t — ts)) for large 
t — ts- The validity of this assumption was partially examined by Sharpe et.al. [1]. Writing 

i?o(t) = Z-^l-AEo-{t-ts) + ^- {AE'.f ■ {t - tsf + 0{{t - tsf)^ , (26) 

they showed in time-ordered perturbation theory that the lattice value of AEq is related to 
the scattering length by the Liischer relation (3) up to corrections of 0{L^^). By a similar 
calculation, one easily shows that the value of AE'q deviates from AEq by terms of 0{L~^). 
These effects occur due to intermediate off-shell two-pion states. 

In the context of our analysis, the momentum cut-off dependence is negligible as discussed 
in Sec. IV. This means that the effects due to the intermediate off-shell two-pion states are 
neghgible. Thus the correction of 0{L~^) for AEq and AE'q is sufficiently small, and the 
time behavior can be regarded as a single exponential function in our simulation. 

To check this point more explicitly, we calculate the scattering length with the energy 
shift obtained with both the linear and the single exponential function in f as was done 
by JLQCD. Results are tabulated in Table VI, which shows that the two set of values are 
consistent within statistical errors, and have no volume dependence. These facts indicate 
that the deviation of JLQCD results between the two fitting functions comes from the 
approximation of the exponential function to the Unear function in t, i.e. the value of 
AEq • {t — ts) ^ / Li^ ■ {t — ts) is not small enough to justify such an approximation due to 
small lattice sizes. 

Another comment concerns the quenching effect on the ratio Ro{t). Bernard and 
Golterman derived the same time behavior (26) using quenched chiral perturbation the- 
ory (qCHPT) [19]. They predicted that the scattering length obtained with quenched ap- 
proximation is divergent in the chiral limit as ao ~ These effects are attributed to 
non-unitarity of the quenched theory. The same results are also obtained by Colangelo and 
Pallante [20]. Divergence in scattering lengths in the chiral limit can also occur if one uses 
chirally non-symmetric lattice fermion action, for example the Wilson fermion action. 

In Fig. 8 we do not observe signs of divergence toward the chiral limit. We consider that 
the effects of quenching and broken chiral symmetry are still too small to affect data at our 
simulation points. 
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The quenching problems can also occur for non-zero momenta, i.e. it is not proven that 
the pion 4-point function Gnmit) behave as a multi-exponential function in i as (15) and 
the diagonalization method can be used. In this work we assume that such effects are small 
in our simulation points as confirmed for the zero momentum case. Investigation of the 
quenching effects for the scattering length and the phase shift by lattice simulations with 
small quark masses is an important future work. 



C. Results of scattering phase shift 

The energy shift AEn = En — En and the phase shift S{pn) at our simulation points are 
tabulated in Tables II, III, and IV. The scattering amplitude A{pn) defined by (25) are also 
included in these tables. 

In Fig. 9 wc plot the amplitude at fixed quark mass as a function of the momentum p^. 
In order to obtain the scattering phase shift for various momenta at the physical pion mass, 
we extrapolate our data with the following fitting assumption : 

p 2 

= Aoo + Aio ■ (ml) + A20 ■ [mlf + Aqi • {f) 

+Ai-(m^)(p2) + Ao2-(py . (27) 

Here A^q corresponds to a^jniT^. In Fig. 9 we omit data plotted with open symbols in the 
fitting. They are for the momentum n = 2 on a 32'^ lattice at niT^jmp = 0.491 and 0.593 
for which a clear plateau in Dn{t) is absent. It should be noted that the constant term Aqq 
vanishes if the effects of quenching and chiral symmetry breaking are negligible. Wc tried 
to fit our data both with and without the assumption Aqq = 0. The results, tabulated in 
Table VII, show that the latter fit yields a value of Aqq which is 1.7a away from zero. The 
other parameters, such as Aiq which are physically more relevant, are consistent between the 
two types of fits, however. From these observations we adopt the value with the assumption 
of Aqq = 0. The fit curves for this fitting are also plotted in Fig. 9. 

Wc present our results for the phase shift 6{p) at the physical pion mass obtained with 
the fitting (27) with the assumption Aqq = in Fig. 10. The filled points arc experimental 
results [21,22]. The values of the phase shift at several momenta are tabulated in Tabic VIII. 
Our results are 30% smaller in magnitude than the experiments. A possible origin of the 
discrepancy is finite lattice spacing effects. As wc saw in Fig. 7 the JLQCD results for scat- 
tering length show a sizable scaling violation. Hence that of the scattering phase shift cannot 
be considered small. Further calculations nearer to the continuum limit or calculations with 
improved actions are desirable to obtain the continuum result of the phase shift. 
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V. CONCLUSIONS 



We have shown in this work that calculations of the scattering length are possible with 
present computing resources. The quenched approximation we employed has theoretical 
issues regarding the chiral extrapolation. We see no problem, either theoretically or compu- 
tationally, in avoiding this problem by going to full QCD calculations, for the simplest case 
of the 1 — 2 two-pion system. The cases of / = and 1=1, which are richer in physics con- 
tent, are much more difficult from the computational point of view. Algorithmic advances 
are presumably needed to evaluate the box and two-loop diagrams with good precision for 
non-zero momenta, which are needed to extract the two-pion energy eigenvalues in these 
channels. 

Another implication of this work is feasibility of a direct calculation of the K tttt decay 
amplitude using the method of Lellouch and Liischer. Diagonalization of the pion 4-point 
function yields the two-pion eigenstate for non-zero relative momenta, which can be used 
as the final state for the K tttt Green function needed in their method. Executing this 
program for the 1 = 2 channel would be an interesting step to take to solve this long-standing 
problem. 
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FIG. 1. Momenta in units of GcV^ used in tliis work for each lattice size. We obtain scattering 
length and phase shift at the filled symbols. Momenta marked by open symbols are used only to 
examine the momentum cut-ofF effects. The broken line shows the upper limit of elastic scattering 
for the smallest m,r in this work, i.e. E = 2 ■ -^/m^T+p^ < 4m7r. 
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FIG. 2. Examples of effective mass of pion propagator at k = 0.1589 {niT^/mp = 0.491) 

on a 32^ lattice. The subscript n refers to the momentum = (27r/L)^ • n. The source is located 
at t = 8. 
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FIG. 3. Examples of the pion 4-point function Gnmit) at k = 0.1589 {m-^jmp = 0.491) on a 
32^ lattice. The two subscripts n and m refer to the sink and source momenta = (27r/L)^ • n and 
k"^ = (27r/L)^ ■ m. The source is located at t = 8. Filled and open symbols indicate positive and 
negative values. In the lower frame, large statistical errors are for Gi2(i), while those of G2i(i) are 
very small. 
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FIG. 4. Ratio Rn{t) and Dn{t) for n = for all quark masses and lattice sizes in this work. 

Quark mass increases from top to bottom, while lattice size increases from left to right. For 
diagonalization of M(i, io), the momentum cut-off is set at iV = 1, and the reference time at 
to = 18. 
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FIG. 5. Ratio Rn{t) and Dn{t) for n = 1 for all quark masses and lattice sizes in this work. 

Quark mass increases from top to bottom, while lattice size increases from left to right. For 
diagonalization of M{t,tQ), the momentum cut-off is set at iV = 1 and 2, and the reference time 
at to = 18. 
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FIG. 6. Ratio Rn{t) and Dn{t) for n = 2 for all quark masses and lattice sizes in this work. 

Quark mass increases from top to bottom, while lattice size increases from left to right. For 
diagonalization of M{t,tQ), the momentum cut-off is set at iV = 2 and 3, and the reference time 
at to = 18. 
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FIG. 8. Comparison of our results on 24^ and 32^ lattices with those of JLQCD on a 16^ lattice 
at 13 = 5.9 [6]. 
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FIG. 9. Scattering amplitude A{p) = taiv5{p)/p ■ E/2 for fixed quark masses. The fit curve is 
also plotted. The open symbols indicate data omitted in the fitting procedure. 
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TABLES 



j 


n = {Nn = 1) 


n = 1 (7V„ = 6) 


n = 2 {Nn = 12) 


1 


-8.913632922 


-1.211335686 


-5.096565798 


2 


16.532315957 


23.243221879 


25.661192388 


3 


8.401923974 


13.059376755 


4.254135936 


A 


6 94^807997 


13 731914368 


14 867^99887 


5 


6.426119102 


11.308518083 


2.283549584 


6 


6.202149045 


13.140942288 


14.148854520 


7 


6.098184125 


11.067054131 


2.051601110 


8 


6.048263469 


13.032596991 


14.031382623 


9 


6.023881707 


11.016034293 


2.011078709 


10 


6.011862830 


13.007939537 


14.007265604 


TABLE L 


Values of the zeta function Zoo(. 


j;n) and = Emez3'^("i^ 


— n) for momenta 



p2 = (27r/L)2 . n. 
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n = (J.iooy 


n 1 FCOQ 

K = U.lOOO 


K = U.lo /4 


K = U.lODO 




(LreV ) 








U. (olyl) 






u.zDUzD(yoj 


0.40896(91) 


11.5405(11; 


V = 24 












Fitting Range 




io — 4U 


1 o An 

18 — 40 


1 O A A 

lo — 44 


1 O A A 

18 — 44 


A 

AEn 


(xltJ GeVj 


98(14) 


97.9(83) 


86.3(53) 


75.5(47) 


Pn -Pn 


(^XlU (oreV j 


oy. 0(^00 j 


en 1 ('/I '^^ 


00.4(^00 j 


OO.y(^oO j 


71 — Tl 


(^XiU J 




1 Q 7"! 


91 fin "i^ 


91 S/'IzL'l 


Pn 


^ X lu V j 




OU. 1 y^o ) 


00.4^^00 J 


oo.y ^^oo j 


"{Pn) 


(deg.) 


— 1 Q1 fS?"! 




— ^ 0'^('95') 


— ^ n7f9fi'l 


A(n ) 




— n 91 4('97'> 


— '^'^5f95'l 


—0 4^8('94'l 


—0 ^'^9('9q'l 


A(fi \ I Tn^ 
■^\Pn)l ''*7r 


l^l/VcreV ) 


— 1 '^4('17'l 


— 1 9Q'^l'Qfi'l 


— 1 11Q('5Q'l 

± . ± ± C/ I tJC/ 1 


—0 97^^51 

U.t/i *JltJ± J 


V = 61 












Fitting Range 




1 O A A 

18 — 44 


1 O A A 

lo — 44 


lo — 44 


18 — 44 


A 77" 


(xlU GeVj 
^'vin— 4 r'Q"\/^2\ 

(^XlU ) 


ol.o(77j 


3o.d(31) 


33.0(21) 


29.3(16) 


Pn -Pn 




1 Q 7/'1 

iy. 1 (i^o) 


91 1 ("1 /I ^ 


91 7^'19^ 


ft — Tl 


(XW ) 


8 7('99'l 


1 7('1 1 
io. / ^lij 


1 A f\f\(Q'^\ 
14.DD^yO j 


1 l^ n9^'89'\ 
10. uzt^oz j 


pi 


(^XiU *jreV ) 


1 9 P.('i^ \ 


iy. 1 1,-LDj 


91 1 (1 A\ 
Zi.H^i4 ) 


91 7^'19"\ 
Zi. ^i-^ j 


0\Pn) 


(deg.) 


U.OUl O-L y 


-L.UZl ±0 I 


1. / Ol -LU J 


— 1 MM 41 


A(n ) 




— n 1 TD^^Q"! 


—0 ^95/'94'l 


—0 43^/'95'l 


—0 ^19f95l 


A(fi \ I ttP' 




— 1 05 ('241 


— 1 250fQ1'l 


— 1 OfiOffiTl 


_n Q^fi('4fii 


V = 48 












r ittmg Range 




18 — 44 


1 O /I /I 

18 — 44 


1 O A 

18 — 4 


1 O /I 

18 — 44 


AEn 


(xlU GeVj 
(xlO-^ GeV2) 


11.8(33) 


in n /o 1 \ 

10.0(21) 


8.6(14) 


7.5(10) 


pI - pi 


4.7(13) 


5.1(11) 


5.51(89) 


5.55(76) 


n — n 


(xlO-3) 


7.4(21) 


7.9(17) 


8.6(14) 


8.7(12) 


pI 


(xlO-^ GeV^) 


4.7(13) 


5.1(11) 


5.51(89) 


5.55(76) 


HPn) 


(deg.) 


-0.68(27) 


-0.74(23) 


-0.84(20) 


-0.85(17) 


MPn) 




-0.217(57) 


-0.294(58) 


-0.399(60) 


-0.464(59) 


A{pn)/ml 


(l/GeV2) 


-1.35(36) 


-1.13(22) 


-0.98(15) 


-0.85(11) 



TABLE II. Results for n = with momentum cut-off N = and to = 18. The scattering 
amphtude A{pn) is defined by A{pn) = taxi. S{pn)/pn ■ En/2. A{pn)lrn^ corresponds to ao/m^. 
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n 1 Q o 
K = U.lOOO 


K = O.io f4 


K = 0.1o6d 




m,-„/mp 






O.D9z(ij 


0. /5z(ij 








U.zDUZD^yUj 


0.40896(91) 


0.o4d8(11 ) 
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Fitting Ranj 
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X iU (jG V j 
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X lU vjrc V J 
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99 ^f 1 4"! 
ZZ. Ol 14 J 


n — n 


(^XiU j 




IZ.Zi 10 j 


Q ^7^ \ 
y.04i^^ / 1 j 
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Cxin~"2 OpV^i 

I /\ XL/ vjrc V J 


9Q 77(7'^') 


98 TfifSQ"! 


98 08('18'l 


97 89^14"! 


"{Pn) 


i^ticg.j 




— 9*^ sf^n"! 


— 18 5('14'l 


— 1 7 Of! 1 "1 

X 1 .UIXX I 
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— n 77(^A^ 


—0 61 1 (7Q^ 


—0 593f4n'l 


—0 594f34'l 


V = 












Fitting Ranj 




io — 00 


io — 40 


1 O /I /I 

18 — 44 


18 — 44 




^^ X lu yjQ V j 


ZZ. O^^iO ) 


1 / .oy (^oi j 


1 c 1 1 (' c:n^ 


lO.U0(oZ ) 


-2 2 


y X iu vjG V J 


IZ.Oo^^oO J 


±1.4Zl oz t 


11. z / \^o0 j 


1 n 88^9^^ 
lU.oo^^ZD J 


Ti — n 


(XiU ) 


a 70(^Q\ 

o. / z^^oy ) 


1 .yZ(^oD J 


( .oZ(^Z4 j 


( .o'iy lo j 


-2 
Pn 


I /N ±\J VJC V J 




1 5 569('59') 


1 5 548('35'l 
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0{Pn) 


l^aeg.j 


— 1 QM 9'l 

-LU. y I -LZ( ) 


— 1 5 9Q('71 
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A(ft 
^\Pn) 








_n f;i Qf 1 K") 

U.OXOl lU J 


U.fJOX I X^ I 


T 7" A o3 

K = 48 












Fitting Ranj 




lo — 44 


1 O A A 

18 — 44 


18 — 44 


18 — 44 




(xlO-3 GeV) 


6.24(81) 


6.08(43) 


5.10(24) 


4.31(22) 


pI-pI 


(xlO-3 GeV^) 


2.96(39) 


3.46(25) 


3.51(16) 


3.37(18) 


n — n 


(xlO-2) 


4.62(60) 


5.40(39) 


5.48(25) 


5.25(27) 


pI 


(xlO-2 GeV^) 


6.705(39) 


6.755(25) 


6.760(16) 


6.746(18) 


5{Pn) 


(deg.) 


-8.8(12) 


-10.35(76) 


-10.50(49) 


-10.06(53) 


A{v„) 




-0.285(38) 


-0.102(29) 


-0.192(23) 


-0.535(28) 



TABLE III. Results for n = 1 with the momentum cut-off N = 1 and to = 18. The scattering 
amplitude A{pn) is defined by A{pn) = tand{pn)/pn • En/2. 
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n 1 c o o 


K = U.io74 


K = 0.15o6 












0. lo2[l) 








U.zDUZD^yUj 


U.4(JoyD(ylj 


0.o4d8(11 j 
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X iU KjG V j 


4U.y(^0D ) 
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Pn Pn 


X lU KjKi V J 


97 S(''?S'l 
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n — n 


(^XiU j 


1 Q "^(07^ 


1 7 /i^ 

i / .U(^14i ) 


1/1 in/'79^ 




Pn 


I /\ XL/ vjrc V J 




31 "^0(90^ 
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30 75fi('7Q'l 


"{Pn) 
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— 18 47fQ3'l 


— 1 7 41 f79'l 


A(n ) 




— n ^7(^(R^ l 


—0 559^46 "1 


—0 507f9fi'l 


—0 590f99'l 


V = 4o 












Fitting Ranj 




io — 00 


18 — 44 


18 — 44 


1 O A A 

lo — 44 




(xlO-3 GeV) 


15.9(10) 


11.40(61) 


9.20(38) 


7.81(33) 


-2 2 
Pn-Pn 


(xlO-3 GeV^) 


8.57(55) 


7.11(38) 


6.75(28) 


6.41(27) 


n — n 


(xlO-2) 


13.37(86) 


11.09(59) 


10.53(44) 


10.00(42) 


-2 
Pn 


(xlO-2 GeV^) 


13.675(55) 


13.529(38) 


13.493(28) 


13.459(27) 


HPn) 


(deg.) 


-17.5(11) 


-14.56(77) 


-13.83(58) 


-13.15(55) 


A{Pn) 




-0.464(31) 


-0.442(24) 


-0.493(21) 


-0.524(22) 



TABLE IV. Results for n = 2 with the momentum cut-ofF N = 2 and to = 18. The scattering 
amphtude A{pn) is defined by A{pn) = taxi. S{pn)/pn ■ -E'n/2. 
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ao/m^ (1/GeV^) qq ■ uIt, 



JLQCD ( LIN ) -2.07(24) -0.0406(47) 

JLQCD ( EXP ) -2.09(35) -0.0410(69) 

Liu et.al. ( Scheme I ) -1.75(38) -0.0342(75) 

Liu et.al. ( Scheme II ) -2.34(46) -0.0459(91) 

CHPT -2.265(51) -0.0444(10) 

TABLE V. Recent results for the scattering length oq in the continuum limit. CHPT refers to 

the prediction of chiral perturbation theory. The error for this case shows theoretical uncertainties. 







K = 0.1589 


K = 0.1583 


K = 0.1574 


K = 0.1566 






0.491(2) 


0.593(1) 


0.692(1) 


0.752(1) 




ml (GeV2) 


0.16113(97) 


0.26026(90) 


0.40896(91) 


0.5468(11) 
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LIN 


-1.23(14) 


-1.194(82) 


-1.042(51) 


-0.917(46) 




EXP 


-1.34(17) 


-1.293(96) 


-1.119(59) 


-0.975(51) 


32^ 


LIN 


-1.02(23) 


-1.207(85) 


-1.029(58) 


-0.912(43) 




EXP 


-1.05(24) 


-1.250(91) 


-1.060(61) 


-0.936(46) 


48^ 


LIN 


-1.34(35) 


-1.12(22) 


-0.97(14) 


-0.84(11) 




EXP 


-1.35(36) 


-1.13(22) 


-0.98(15) 


-0.85(11) 



TABLE VI. Our results for the scattering length ao/mT^ (1/GeV^) calculated from the energy 
shift obtained by the liner fitting (LIN) and the exponential fitting (EXP) of Ro{t) in t. 



Chiral No-Chiral 



^00 




-0.069(41) 


^10 (l/GeV2) 


-1.389(84) 


-1.01(24) 


A20 (l/GeV^) 


0.79(18) 


0.33(33) 


^01 (l/GeV2) 


-2.07(20) 


-2.00(20) 


^11 (l/GeV^) 


3.22(47) 


3.09(48) 


A02 (1/GeV^) 


1.27(53) 


1.23(53) 




0.863 


0.782 



TABLE VII. Results of fitting of the scattering amplitude with the assumption (Chiral), and 
without the assumption Aqq = (No-Chiral). 
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(GeV^ 

0.020 

0.070 

0.140 

0.230 

0.340 



(GeV) 



0.40 
0.60 
0.80 
1.00 
1.20 



S{p) (deg.) 
-2.71(12) 
-8.09(59) 
-14.8(12) 
-22.0(20) 
-28.6(31) 



TABLE VIII. Our results for the scattering phase shift at several momenta at the physical pion 
mass. 
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